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ABSTRACT 

Pressure and velocity data were collected in a full scale jet 
engine test cell in order to validate the predictive accuracy of a 
two-dimensional and axisymmetric primitive variable computer code. It 
was found that the model reasonably predicted the velocity profiles in 
the augmentor tube. Inaccuracy increased at higher engine thrust settings 
at positions far downstream in the augmentor tube. Predicted pressure 
profiles were reasonable but the magnitudes were in considerable error 
at high flow rates. 
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I. INTRODUCTION 



Jet engine test cells are in wide use today in military and com- 
mercial aviation maintenance programs. Their purpose is to provide 
for testing of jet engines under conditions that approximate the in- 
stalled environment, thus enabling the responsible authority to veri- 
fy the engine's capabilities throughout its operational envelope 
prior to returning it to service. However, due to exhaust pollution 
problems inherent to their operation, the use of these test facilities 
has recently encountered a considerable amount of criticism and even 
judicial action in the State of California. 

The necessity for testing an engine after overhaul is apparent. 
However, the pollution problems that result when testing today's high 
power, high mass flow engines in present test facilities must also be 
considered. 

"The typical test cell incorporates an inlet, a horizontal test 
section and vertical exhaust stack. The engine to be tested is 
normally mounted near the center of the cell to allow the development 
of a nearly uniform engine inlet velocity profile. The engine ex- 
hausts into an augmentor tube which entrains additional air for ex- 
haust gas cooling and dilution. The quantity of this secondary air 
is crucial to proper engine testing and test cell performance" 

[Ref. 1]. The ratio of entrained air mass flow rate to engine mass 



10 



flow rate is known as "augmentation ratio". An augmentation ratio that 
is too high could result in inaccurate engine performance measurements 
and in exceeding the test cell structural limits due to pressure 
gradients and excessive cell pressure reduction. Insufficient secondary 
air could allow engine exhaust gas recirculation to the engine inlet and 
cause excessive opacity of visible emittants. It is therefore necessary 
to develop mathematical models which can be used to predict the flow 
fields within turbojet test cells. These models could then be used as 
a cost effective means for determination of optimum test cell designs 
for meeting the military engine testing needs while minimizing noise and 
chemical pollution effects on the local environment. 

With this in mind, considerable experimental and computer modelling 
work has been done at the Naval Postgraduate School to further the 
understanding of the flow field and detailed operating characteristics 
encountered in typical test cells. This knowledge is currently needed 
to facilitate the necessary modifications to present test cells and for 
the design of new cells to accommodate the high technology engines of 
the future. 

The two-dimensional modelling of the flow within turbojet test 
cells was initiated at NPS by Hayes and Netzer [Ref. 2]. They adapted 
the two-dimensional, ip-oj model of Spalding, et al [Refs. 3, 4] to the 
test cell geometry and operating conditions. This work was continued 
by Speakman, Hayes and Netzer [Ref. 5]. The model made use of the 
fundamental equations of conservation of mass, momentum and energy. 

The dependent variables were stream functions (tp) , vorticity (co). 
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turbulence kinetic energy (K), the turbulence dissipation rate (e), and 
the temperature (T). As a consequence of the ip-oj formation, pressure 
(and to a large extent velocity) is removed from the governing equa- 
tions. These second order elliptic partial differential equations 
were reduced to "finite-difference," non-linear algebraic equations and 
solved simultaneously with an iterative procedure. Results of this 
computer simulation appeared to be quite accurate at low subsonic en- 
gine exhaust Mach numbers (M < 0.6), but pressure predictions were un- 
reliable at high subsonic or sonic exhaust velocities. These compari- 
sons between predictions and experiment were made by Walters and Netzer 
[Ref. 6] using a one eighth scale model of an NAS Alameda test cell. 
Other problems related to the ip-u formulation were: 1) restriction to 
constant density flows, or to flows in which density varied only with 
temperature; 2) boundary conditions were difficult to specify; 3) con- 
siderable difficulty was experienced in obtaining converged solutions, 
especially for non-uniformly spaced grids and high flow rates, and 4) 
the ip-co model is not easily extended to three dimensional flows. 

To help alleviate the above modelling difficulties Stevenson and 
Netzer {Ref. ^utilized a primitive variable ( pressure- velocity) model . 
This computer model (adapted from the CHAMPION 2/E/FIX computer program 
developed by Pun and Spalding ^Ref. 7^ extended the modelling capabili- 
ties to more complex geometries and to high subsonic and sonic engine 
exhaust velocities. It also is readily applicable to variable density 
flows. The earlier results of the i/j-uj model and those of the u-v-p 
model were compared to the empirical data obtained with the subscale 
test cell {jtef. 6*]. Good agreement between the two methods was found 
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for low subsonic engine exhaust conditions. The primitive variable 
computer model produced reasonable results for the flow field downstream 
of the engine exhaust up to the sonic exhaust condition. 

The primary objective of this investigation was to adapt Spalding's 
CHAMPION 2/E/FIX primitive variable computer model (as modified by 
Stevenson) to the geometry and flow conditions of an actual operational 
full scale turbojet test cell and to test that model with empirical data 
from the same cell. The turbojet test cell chosen was one located at 
NAS Alameda, operated by NARF Alameda. The engine being operated at 
the time of data acquisition was the General Electric TF34-GE-100, cur- 
rently used on the USAF A-10 attack aircraft. The major computer model 
modifications required were: 1) to increase the geometry to the full 

scale, and 2) to change the uniform engine exhaust velocity profile to 
a non-uniform profile which adequately represented the primary and 
secondary flows from the TF-34 turbofan engine. 
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II. EXPERIMENTAL INVESTIGATION 



A. INTRODUCTION 

The turbojet test cell chosen for this investigation was test cell 

No. 13 located in Building 372 at NAS Alameda. It is operated by NARF 

Alameda as the prime means of evaluating the USAF TF34-GE-100 turbofan 
engine following required normal periodic maintenance. 

The TF34-GE-100 is a turbofan engine in current use on the USAF 
twin engine A- 1 0 attack airplane. The engine is 2.54m long and produces 
a maximum rated thrust of 9065 lbf JjRef . 8j. The bypass ratio is 6.22:1 
with maximum rated airflow at 151.0 kg/sec. The exhaust nozzle is a 
fixed area, convergent type with the cold stream being 0.474 sq. m. in 
area and the hot stream 0.148 sq. m. in area. While operating at maxi- 
mum power the core exhaust gas is estimated to have a velocity of 305 m/ 

sec and a temperature of 393°C jjRef. 8^. 

To provide the data for model validation it was decided to measure 
velocity profiles within the test cell and at two locations within the 
augmentor tube as well as the axial pressure distribution throughout 
the test cell. To most accurately determine velocity profiles, mul- 
tiple pitot-static tubes with associated stagnation temperature measure- 
ments would have been required. However, due to cost and time restraints 
involved with the manufacturing and installation of such equipment, a 
compromise approach to this problem was taken. Two stagnation tube 
rakes, each fitted with two thermocouples, were installed in the aug- 
mentor tube. Also, thirty two static pressure ports were installed in 
the augmentor tube wall. Velocity measurements were also made within 
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the test cell using hot wire velocimeters . Details of the data acquisi- 
tion as well as the method of data reduction are presented below. 

B. TEST CELL 

The turbojet test cell utilized in this investigation is typical of 
test cells currently used in military and commercial aviation maintenance 
programs. It was a concrete blockhouse type structure 43m long, 9m high, 
and 9m wide (Fig. 1). It had a horizontal inlet, which was retangular in 
cross section, and exhausted into a vertical stack. Not shown in Figure 
1 are the noise suppression devices which were located over the cross 
sections of both the test cell inlet and back wall. Located 15m from 
the cell inlet was the test platform which stood 3m above the floor of 
the cell. The jet engine test stand and the first section of the aug- 
mentor tube were located on this platform. The TF-34 was mounted on the 
test stand with an attached, fixed geometry bellmouth type inlet, along 
with necessary start-up and test equipment to monitor engine performance 
in the cell control room. This test stand was movable, so that the 
distance between the exit plane of the engine and entrance plane of the 
augmentor tube could be varied. During this investigation this distance 
was fixed at 1 .7m. 

The first section of the augmentor tube was located on the test 
platform. As shown in Figure 1 it too was movable. It was a steel 
pipe with 9.5mm thick walls and was approximately 1.8m in diameter and 
1.6m long. The next section was also a constant diameter pipe, 5.2m 
long. During this investigation, the pipe protruded only 5m beyond 
the first section. The radial gap between these two sections was, on 
the average, 16mm. The next section was 3.5m long and diverged in 
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diameter from approximately 1.8m to 2.7m. At the downstream end of the 
latter section was a water-injection ring. Following a 0.31m gap was 
the next section: a constant 3.05m diameter pipe, 13.1m long. In 

this section were three staggered rows of four noise suppression de- 
vices. Each unit was a 3.2m long cylinder, 0.6m in diameter. 

C. INSTRUMENTATION AND DATA ACQUISITION 

Two pitot tube rakes were utilized in this investigation. Each 
rake was constructed from 102mm wide, 9.5mm thick, 302 stainless steel. 
The cross-like configuration used is depicted in Figs. 2 and 3. The 
first rake was located 127mm from the entrance plane of the augmentor 
tube. The tips of the pitot tubes being 25.4mm from the entrance. 

The vertical member of this rake was 1.67m high and had seven pitot 
tubes built into it. Each horizontal member was 0.822m long with three 
pitot tubes inserted. Horizontal and vertical members were attached as 
shown in Fig. 4. To hold the rakes in position, 50.8mm high tabs were 
welded into the augmentor tube at the 12, 3, 6, and 9 o'clock positions. 
The rakes were then attached to the tabs as shown in Fig. 5. The down- 
stream rake was constructed exactly the same as the first except that 
the vertical length was 1.70m and each horizontal member was 0.833m 
long. This second rake was installed 5.05m downstream from the aug- 
mentor entrance plane (127mm upstream from the start of the diverging 
section of the tube). The 6.35mm diameter stainless steel tubing used 
for the pitot tubes had a wall thickness of 1.57mm. Each tube pro- 
truded 101.6mm forward from the rake. The tubes were silver-soldered 
on the back side of the rake and routed from the augmentor tube at the 
aft end of each attachment tab through holes drilled in the augmentor. 
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Tygon tubing (3/16 in. ID x 1/16 in. wall) was then routed from pitot 
tube to a Scanivalve located in a standing cabinet next to the jet- 
engine test stand. The output of the Scanivalve was routed through the 
floor of the cabinet and underneath the test platform to the recording 
apparatus in the control room (Fig. 6). In addition to the pitot tubes, 
thermocouples were mounted on each rake as shown in Fig. 2. The thermo- 
couple wire was routed from the augmentor with the pitot tubing, then to 
ice bath references. The thermocouple data were obtained from a digital 
voltmeter. 

Static pressure ports were located at the 2, 5, 7, and 10 o'clock 
positions at axial locations 0.152m, 0.305m, 0.610m, 1.212m, 3.96m, 
5.18m, 6.10m, and 7.32m downstream of the augmentor inlet. Circumferen- 
tial ports were used to determine whether or not the flow was axi sym- 
metric. Tygon tubing was attached to Swagelock fittings at each port 
and routed to a second Scanivalve in the standing cabinet. 

In order to determine the velocity of the secondary (test cell) air, 
velocimeters were located 0.6m and 1.8m above the TF34 engine. The 
velocimeters used were Datametrics Air Flow Meters, commonly used to 
measure the velocity in air conditioning ducts. Prior to their use, 
calibration was performed at the NPS using a small wind tunnel. The 
resulting- calibration curve is shown in Fig. 7. The output of these 
velocimeters were routed underneath the jet engine test stand to the 
standing cabinet and then to the control room as described above. In 
the control room the output signals were connected to a digital volt- 
meter. 

The instrumentation apparatus set up in the control room (Fig. 6) 
consisted of a locally manufactured Scanivalve controller, two Doric 
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digial voltmeters and two Honeywell "Electronic 196", 2 pen strip chart 
recorders . 

D. EXPERIMENTAL PROCEDURE 

Prior to actual data acquisition the following steps were performed: 

1) The 30.5 cm radial gap between the second and third augmentor 
sections was blocked off with 1.9 cm thick plywood. This was 
done to prevent any secondary cell air from entering the aug- 
mentor downstream of the data acquisition equipment. It should 
be noted that an airtight seal was not attempted. 

2) All static and dynamic pressure ports were pressure checked to 
insure proper connections and identification. 

3) All thermocouples were checked for proper rake attachment, 
wiring and identification. 

The actual data acquisition process began before anticipated due to 
inadequate thermocouple installation. After initial TF34 start up and 
a brief run time (during which the recording instruments were calibrated) 
the test cell operators discovered an oil leak in the TF-34 and conse- 
quently shut it down for repair. This provided an opportunity to 
inspect the test equipment. At this time it was discovered that the 
thermocouples were not suitably constructed to withstand a complete 
series of tests. Therefore, it was decided to record termperature data 
during the initial TF-34 bearing breakin tests at various power settings 
for as long as the thermocouples remained in place. Reliable tempera- 
ture data were only obtained for maximum, 90%, and idle power settings. 
The data appears in Table I. 
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POWER 

SETTING 



100 % 



90% 



IDLE 



TABLE I: 



CORE FAN THRUST RAKE 
RPM RPM LBS 

FORWARD 

17228 6667 8611 

AFT 

FORWARD 

16120 5803 5967 

AFT 

FORWARD 

11350 1685 382 

AFT 



THERMOCOUPLE 


TEMPERATURE 

Cc) 


CENTER 


438 


LOWER 


39 


CENTER 


127 


LOWER 


76 


CENTER 


461 


LOWER 


32 


CENTER 


156 


LOWER 


62 


CENTER 


498 


LOWER 


18 


CENTER 


105 


LOWER 


42 



Thermocouple Data Obtained During Initial Engine Testing 
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After the brief bearing breakin tests the TF-34 was put through a 
series of five minute duration tests at a variety of power settings. 
When a particular power setting was reached two minutes were allowed 
to insure steady-state operation and then the following data were re- 
corded: 

1) From the test cell instrumentation: 

a) ambient temperature 

b) ambient pressure 

c) cell pressure 

d) TF-34 core RPM 

e) TF-34 fan RPM 

f) TF-34 thrust 

2) From NPS installed equipment and instrumentation: 

a) static pressures along the augmentor wall 

b) pitot tube pressures at each rake 

c) temperatures at each rake 

Each 48 channel Scanivalve was set up so that port number one was 
open to the cell pressure. A fixed pressure (set using a mercury 
manometer) was input to port number two in order to provide a refer- 
ence calibration each time a cycle was begun. During each run one 
Scanivalve was cycled using the Scanivalve controller, data recorded 
on the strip chart, connections switched, and the same procedure 
followed for the remaining Scanivalve. The pitot rakes were labeled 
as shown in Fig. 8, and static ports as in Fig. 9. The recorded data 
are presented in Tables II through IV, the letter "A" referring to the 
forward rake, "B" to the aft rake, and letters "C through J" to the 
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static port stations from front to rear. The data obtained from the 
intermediate power setting run was not utilized due to lack of reliable 
temperature information at that thrust level. 

There were three specific problems encountered during the data ac- 
quisition process: 

1) Failure of thermocouple attachments. This resulted in lack of 
temperature data during each run that pressure data were ob- 
tained. However, this was overcome by correlating (by thrust 
setting and fan and core speeds) the temperature data recorded 
at the start of the TF-34 test series with the pressure data 
obtained later. 

2) During the entire test series a loud resonant frequency was 
heard. This initially forced engine shut down to allow for 
cell inspection. No damage was discovered and all instrumen- 
tation appeared normal so testing was continued as planned. 
However, upon removal of test equipment it was discovered that 
two of the welds holding the aft rake in place had broken. It 
should also be noted that at the same time it was discovered 
that two pitot tubes had detached from the same rake. The 
latter apparently occurred during the last test since no 
unusual change occurred in the recorded data. 

3) Due to electrical connection difficulties the velocimeters were 
not installed during the engine testing periods when other data 
were obtained. When this problem was corrected (after the pres- 
sure data were obtained) it was discovered that the wooden blocks 
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RUN NO.: 1 



POWER SETTING: 100% 



TIME: 1305 



PITOT RAKE (P-P CELL ) STATIC PORT (P-P CELL ) 



SCANIVALVE 

CHANNEL 


PORT I.D. 


(N/m 2 xl0 -3 ) 


SCANIVALVE 

CHANNEL 


PORT I.D. 


(N/m 2 xl0 -3 ) 


1 


Zero 




1 


Zero 




2 


REF PRES 




2 


REF PRES 




3 


A1 


34.30 


3 


Cl 


-14.196 


4 


A2 


39.93 


4 


C2 


-13.776 


5 


A3 


3.56 


5 


C3 


-13.803 


6 


A4 


-0.63 


6 


C4 


-14.266 


7 


A5 


41.11 


7 


D1 


-13.218 


8 


A6 


5.97 


8 


D2 


-12.797 


9 


A7 


-0.48 


9 


D3 


-12.956 


10 


A8 


39.58 


10 


D4 


-13.190 


11 


A9 


3.78 


11 


El 


-9.694 


12 


A10 


-0.84 


12 


E2 


-9.370 


13 


All 


40.19 


13 


E3 


-9.432 


14 


A12 


2.68 


14 


E4 


-9.860 


15 


A13 


-0.48 


15 


FI 


-8.226 


16 


B1 


26.84 


16 


F2 


-8.226 


17 


B2 


19.81 


17 


F3 


-8.226 


18 


B3 


8.39 


18 


F4 


-8.226 


19 


B4 


3.30 


19 


G1 


-6.433 


20 


B5 


23.24 


20 


G2 


-6.661 


21 


B6 


10.36 


21 


G3 


-6.530 


22 


B7 


3.12 


22 


G4 


-6.433 


23 


B8 


18.71 


23 


HI 


-6.268 


24 


B9 


7.51 


24 


H2 


-6.171 


25 


BIO 


1 .28 


25 


H3 


-5.875 


26 


Bll 


15.59 


26 


H4 


-6.33 


27 


B12 


5.76 


27 


11 


-4.144 


28 


B13 


0.04 


28 


12 


-4.013 








29 


13 


-3.820 








30 


14 


-4.047 








31 


J1 


-1.731 








32 


J2 


-1 .958 








33 


J3 


-2.027 








34 


J4 


-1 .827 



FAN 

RPM: 6701 


CORE 

17360 






THRUST: 8685 

(LBS) 

CELL 
TEMP: 17.8 

(°c) 


RAKE A 


RAKE B 


CELL AMBIENT 

PRESSURE: 1 .01 29x10^ 1.0141xl0 5 

(N/m 2 ) 



TABLE II: EXPERIMENTAL DATA, 100% POWER 
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RUN NO: 3 

PITOT RAKE (P-P CELL 


POWER : 

) 


SETTING: 90% 

STATIC 


PORT (p- 


TIME: 1327 

P CELL> 


SCANIVALVE 

CHANNEL 


PORT I.D. 


(N/m 2 xlO 


-3x SCANIVALVE 
> CHANNEL 


PORT I.D. 


(N/m 2 xlO -3 ) 


1 


Zero 




1 


Zero 




2 


REF PRES 




2 


REF PRES 




3 


A1 


22.64 


3 


Cl 


-9.29 


4 


A2 


30.86 


4 


C2 


-10.08 


5 


A3 


0.49 


5 


C3 


-9.56 


6 


A4 


-0.39 


6 


C4 


-9.29 


7 


A5 


29.04 


7 


D1 


-9.23 


8 


A6 


4.18 


8 


D2 


-9.23 


9 


A7 


-0.19 


9 


D3 


-9.33 


10 


A8 


30.60 


10 


D4 


-9.45 


11 


A9 


2.58 


11 


El 


-7.46 


12 


A10 


-0.29 


12 


E2 


-6.70 


13 


All 


31 .71 


13 


E3 


-7.07 


14 


A12 


3.65 


14 


E4 


-7.55 


15 


A13 


-0.26 


15 


FI 


-5.71 


16 


B1 


20.64 


16 


F2 


-5.62 


17 


B2 


16.86 


17 


F3 


-5.53 


18 


B3 


8.18 


18 


F4 


-5.59 


19 


B4 


1.60 


19 


G1 


-4.14 


20 


B5 


17.30 


20 


G2 


-4.41 


21 


B6 


8.41 


21 


G3 


-4.41 


22 


B7 


3.16 


22 


G4 


-4.26 


23 


B8 


13.96 


23 


HI 


-4.08 


24 


B9 


5.83 


24 


H2 


-3.87 


25 


BIO 


1.38 


25 


H3 


-3.75 


26 


Bll 


12.72 


26 


H4 


-4.12 


27 


B12 


3.96 


27 


11 


-2.76 


28 


B13 


0.98 


28 


12 


-2.79 








29 


13 


-2.45 








30 


14 


-2.61 








31 


J1 


-1 .10 








32 


J2 


-1.40 








33 


J3 


-1.40 








34 


J4 


-1.34 



FAN 


CORE 






RPM: 5933 


16403 






THRUST: 6304 






(LBS) 








CELL 


RAKE A RAKE B 


CELL , 


AMBIENT 


TEMP: 18.3 

(°C) 




PRESSURE: 1 .01 36x1 0 D 

(N/nr ) 


1 .0144x10 



TABLE III: EXPERIMENTAL DATA, 90% POWER 
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RUN NO: 4 




POWER SETTING 


: IDLE 




TIME: 1337 


PITOT 


RAKE (P“Pq£[_[_ 


) 


STATIC PORT (p- 


P CELL^ 


SCANIVALVE 

CHANNEL 


PORT I.D. 


(N/m 2 xl0 -3 ) 


SCANIVALVE 

CHANNEL 


PORT I.D. 


(N/m 2 xl0~ 3 ) 


1 


Zero 




1 


Zero 




2 


REF PRES 




2 


REF PRES 




3 


A1 


1.44 


3 


Cl 


-0.64 


4 


A2 


2.56 


4 


C2 


-0.77 


5 


A3 


0.45 


5 


C3 


-0.74 


6 


A4 


0.27 


6 


C4 


-0.71 


7 


A5 


2.52 


7 


D1 


-0.71 


8 


A6 


0.41 


8 


D2 


-0.74 


9 


A7 


0.27 


9 


D3 


-0.74 


10 


A8 


2.43 


10 


D4 


-0.71 


11 


A9 


0.27 


11 


El 


-0.58 


12 


A10 


0.14 


12 


E2 


-0.49 


13 


All 


2.48 


13 


E3 


-0.52 


14 


A12 


0.23 


14 


E4 


-0.55 


15 


A13 


0.14 


15 


FI 


-0.40 


16 


B1 


1 .53 


16 


F2 


-0.43 


17 


B2 


1 .48 


17 


F3 


-0.40 


18 


B3 


0.54 


18 


F4 


-0.43 


19 


B4 


0.50 


19 


G1 


-0.28 


20 


B5 


0.27 


20 


G2 


-0.28 


21 


B6 


1.35 


21 


G3 


-0.25 


22 


B7 


1.35 


22 


64 


-0.28 


23 


B8 


0.36 


23 


HI 


-0.28 


24 


B9 


1.12 


24 


H2 


-0.18 


25 


BIO 


0.54 


25 


H3 


-0.18 


26 


Bll 


0.32 


26 


H4 


-0.15 


27 


B12 


1.12 


27 


11 


-0.03 


28 


B13 


0.41 


28 


12 


-0.03 








29 


13 


0 








30 


14 


0 








31 


J1 


-0.03 








32 


J2 


0 








33 


J3 


0 








34 


J4 


0 



FAN CORE 
RPM: 1801 11404 



THRUST: 410 

(LBS) 

CELL RAKE A RAKE B 
TEMP: 18.3 

(°C) 



CELL , AMBIENT , 

PRESSURE: 1.0144x10® 1. 0144x10® 

(N/m 2 ) 



TABLE IV: EXPERIMENTAL DATA, IDLE POWER 
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at the intersection of the second and third augmentor sections 
had partially detached. This caused the velocity of the cell 
air to probably be somewhat higher during the tests using the 
velocimeters than it was during the pressure/acquisition runs 
when the blocks were in place. The velocimeter data are shown 
in Table V. These data were used only to obtain the general 
shape of the velocity profile above the engine. 



E. DATA REDUCTION 

The (p-p cell ) data that appear in Tables II through IV came directly 
from the strip chart recordings, or, in the case where the reading was 
too small, from the digital voltmeter. The calibration scale on the 
strip chart (psi/division) was obtained from a 1.05 in. Hg reference in- 
put to Channel No. 2 of both Scanivalves. Using this information, 
axial pressure profiles along the four circumferential positions were 
plotted and an average pressure profile was then estimated for each 
power setting (Figs. 10, 11, and 12). These profiles were then used to 
obtain the static pressure at the aft rake location (it was assumed 
that pressure was uniform across the augmentor at the aft rake). 

The measured stagnation pressure at each pitot tube together with 
the static pressure at the rake were used to obtain a Mach number pro- 
file using the expression 




X 



Y 



25 



Velocity RPM Thrust 

(m/sec) (lbs) 



Upper 

Velocimeter 


Lower 

Velocimeter 


Fan 


Core 




0.51 


1.14 


1789 


11265 


520 


2.03 


3.05 


6780 


1 7285 


9072 



TABLE V: VELOCIMETER DATA 
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with y assumed equal to 1.4. The thermocouple readings were plotted and 
used to estimate an assumed axi symmetric stagnation temperature profile 
for the aft rake at each power setting (Fig. 13). The estimated stagna- 
tion temperature profile and local Mach number were then used to deter- 
mine a static temperature for each pitot tube location. Using this in- 
formation and the properties of air at standard temperature and pressure, 
accoustic speed and velocity were found for each pitot tube position. 
Static temperature for the aft rake is plotted as a function of radius 
in Fig. 14, velocity profiles are shown in Figs. 15, 16, and 17. 

The static pressure of the front rake was estimated by extrapolating 
the average pressure profile to the augmentor inlet position. This was 
required due to the fact that the first set of static pressure ports 
was installed 15.2 cm aft of the augmentor entrance plane. This static 
pressure was assumed to be uniform across the augmentor for each power 
setting. Mach number profiles were then obtained as above. At the 
forward rake y was not assumed to be uniform across the augmentor. In- 
stead, y was taken to be 1.33 at the centerline position and at the four 
pitot tubes located circumferentially 25.4 cm from the center. At the 
remaining positions y was assumed to be 1.4. This property change was 
made since heat induced discoloration on the forward rake was observed 
to occur from the core jet to a radius of approximately 30 cm. Stagna- 
tion temperature profiles were estimated using the measured temperatures 
at the two rake locations and are shown in Fig. 18. Using this stagna- 
tion temperature distribution with the Mach number profiles, accoustic 
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speed and velocity were determined as above. The resulting velocity 
profiles for each power setting are shown in Figs. 19 through 21. 
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III. COMPUTER MODEL 



A. INTRODUCTION 

The computer program used to model the turbojet test cell of this 
investigation was the CHAMP ION/2/ E/F IX computer program developed by 



primitive variable (pressure-velocity) model for two-dimensional, turbu- 
lent recirculating flows. The equations are cast into finite difference 
form and solved using a line-by-line iteration method. A detailed dis- 
cussion of the model can be found in Refs. 1 and 7. 

B. PHYSICAL MODEL 

The major modifications made to the computer program as utilized by 
Stevenson ^Ref. ij were to change the geometry of the test cell and to 
incorporate the engine exhaust velocity profile of the TF-34 turbofan 
engi ne. 

The turbojet test cell, TF-34 engine and augmentor dimensions used in 
the computer model are shown in Fig. 22. There were several differences 
between the configuration used in the computer model and the actual test 
cell as shown in Fig. 22: 

1. The computer model was axi symmetric and therefore required the 
engine to be located on the axis of symmetry of the test cell. 

In the actual cell the engine is approximately 0.61m closer to 



the deck than to the overhead. Thus, it can be expected that the 
predicted velocity distribution of the secondary cell air would 
be somewhat different than the actual distribution. In addition. 



Pun and Spalding 




Stevenson 




This is a 
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the test stand and supporting deck probably distorted the velocity 
profile at the augmentor inlet. 

2. The exhaust of the TF-34 turbofan engine was modeled as two con- 
centric cylinders (gas generator core and fan exhausts). In 
reality, each exit incorporates a converging nozzle. For simpli- 
city in the two-dimensional model, the actual exit diameter of 
each duct was utilized over the entire engine length (Figs. 23 
and 24). 

3. In the computer model the aft cell wall was taken to coincide 
with the augmentor inlet. This wall was actually approximately 
11m farther back. The effect of this modification was to elimi- 
nate the reci rculation zone above the augmentor tube. However, 
the effects of this low velocity reci rculation zone above the 
augmentor flow field have been shown to be insignificant [Ref. 5], 

4. The computer model utilized a constant diameter augmentor tube 
20.5m long instead of the varying area tube construction shown 
in Fig. 1. The second rake was not positioned farther down- 
stream to avoid complications arising from the divergent 
character of the second augmentor section and the flow distur- 
bance expected to result from the water injection ring and 
noise suppression devices. The augmentor tube length was 
modelled as 20.5m versus the actual 5.05m so that the boundary 
conditions at the exit plane would not influence the calcula- 
tions at the position of the second rake. 

5. The test cell was square in cross-section. This was modelled as 
a concentric cylinder with the same cross-sectional area. 



30 



The grid spacing for the finite difference numerical solution of this 
problem was crucial. A 30 by 30 grid was used originally without obtain- 
ing convergence. Consequently the grid was expanded to a 40 by 30 to 
better accommodate the much greater length of the axial direction. The 



not increase by more than a factor of approximately 1.5 was adhered to. 

C. DATA INPUT 

Using the empirical velocity profiles from Figs. 15 through 17, an 
average velocity for each radial position at the aft rake was determined. 
These average velocities and the centerline velocity were then plotted 
(Fig. 25) and graphically integrated to obtain a one-dimensional bulk 
velocity at the aft rake for each power setting. An identical procedure 
was followed to obtain an average static temperature (Fig. 14) at the aft 
rake. These results were then used (together with the standard properties 
of air and known augmentor tube cross sectional area) to calculate the 
total mass flow rate through the augmentor tube. 

Using the fan and core RPM, thrust, and engine inlet temperatures 
that appear in Tables II through IV, fan and core mass flow rates for 
the TF-34 were obtained from the General Electric Corp., Lynn, Mass. 

It should be noted that the flow rates provided by General Electric 
were for "an ideal" engine and may not exactly apply to the specific 
engine used in this investigation. 

With the above flow rates it was then possible to calculate the 
mass flow rate through the test cell. Then, using the temperature and 
pressure in the test cell (and the known cell cross sectional area) the 
average one-dimensional velocity of the cell air was determined. 



recommendation by Gosman 




should 
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In this investigation no instrumentation was available for measuring 
static pressures and temperatures at the TF-34 core and fan tailpipe ex- 
hausts. Turbine exit temperature as a function of thrust is presented in 
Ref. 8. 

This temperature was used as the tailpipe stagnation temperature in 
order to calculate nozzle exhaust temperature. For the idle power setting 
the Mach number at the core exit was assumed to be approximately 0.5 and 
y to be 1.33. Using a turbine exhaust temperature of 1560°R resulted in 
a nozzle discharge static temperature of 1498°R. For the 90% and full 
power cases the nozzle flow was choked and the turbine exhaust tempera- 
tures were 1694°R and 1922°R respectively. This resulted in exhaust 
temperatures of 1453°R and 1649°R respectively. Fan exhaust temperatures 
were available from the temperature profiles of Ref. 8. Actual cross 
sectional areas of both exhaust nozzles were measured on the engine. 

Static pressure at the exhaust nozzles was taken to be the static pres- 
sure of the cell in all cases. This static pressure estimate is 
accurate for the fan exhaust in all cases and the core exhaust at the 
idle power setting. However, the static pressure of the core exhaust 
for the maximum and 90% power settings was actually the critical pres- 
sure. 

Using the above information the one-dimensional continuity equation 
was used to calculate the bulk velocities for core and fan exhausts at 
all power settings. The computer input data as well as the CPU time 
required to obtain a converged solution for each power setting are shown 
in Table VI. 
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POWER SETTING 




TEMPERATURE 

(°K) 


VELOCITY 

(m/sec) 


CELL-PRESSURE 

(N/nrxl0~ 3 ) 


CPU TIME 

(IBM 360/67) 
(min) 




CELL 


291 


0.825 


101.5 


75 


IDLE 


FAN 


344 


56.04 








CORE 


832 


59.63 








CELL 


291 


2.75 


101.37 


75 


90% 


FAN 


305 


174.2 








CORE 


807 


258.8 








CELL 


290 


3.52 


101.30 


150 


100% 


FAN 


306 


202.9 








CORE 


916 


354.1 







TABLE VI. COMPUTER INPUT DATA 
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IV. DISCUSSION OF RESULTS 



A. INTRODUCTION 

The results of this investigation are presented in the following 
manner: 

1) General characteristics and accuracy of the experimental data 
and sensitivity of the calculated experimental velocity pro- 
files to assumptions made for the unknown/unmeasured experi- 
mental parameters. 

2) A comparison of experimental data to computer model predictions 
and the sensitivity of the predictions to uncertainties in the 
data. 

B. EXPERIMENTAL DATA 

Experimental data are presented in Tables II through IV. The re- 
sulting calculated velocity profiles for the fore and aft rakes at each 
power setting are presented in Figs. 19 through 21 and Figs. 15 through 
17, respectively. Figs. 10 through 12 show the static pressure profiles 
along the augmentor wall for each power setting. 

The pressure profiles show only a small variation with circumferen- 
tial location. In all cases the minimum pressure occurred very close 
to the inlet plane of the augmentor, whereas the maximum pressure in 
the instrumented portion of the augmentor at no time exceeded atmos- 
pheric. As expected, the minimum pressure at the augmentor inlet de- 

-3 2 

creased as the engine flow rate increased, going from -0.8x10 N/m 

-3 2 

at the idle power setting to -14.0x10 N/m at the maximum power condi- 
tion. 
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The velocity profiles at the aft rake indicate that the flow was not 
axisymmetric in the augmentor tube. The most asymmetric condition oc- 
curred at the idle power setting. The idle profile was also much more 
fully developed than the others, indicating that the flow at idle was 
the only one that entered the diverging section of the augmentor tube in 
a reasonably well mixed condition. There also appeared to be a general 
trend toward higher mass flow in the quadrant defined by the 12 and 3 
o'clock positions. 

In all cases the velocity profiles at the forward rake showed the 
expected high velocity near the center with the lower fan and augmentor 
air velocities toward the wall. It should also be noted that the spread 
in data at the forward rake was much less than that at the aft rake. 

In all cases the engine-exit-to-augmentor-inlet spacing was 1.7m. At 
the idle power setting the augmentation ratio was 2.39:1. This ratio 
dropped to 2.22:1 for the 90% power case and rose to 2.42:1 for the maxi- 
mum power setting. 

In summary it can be said that the static pressure remained quite 
uniform across the tube, that the flow at the aft rake was not axisymmet- 
ric, and that at the front pitot rake the flow was quite symmetrical. 

For the 90% power setting the stagnation pressure data at the aft rake 
was estimated to have an uncertainty of 0.17% which, in turn, resulted in 
a 0.25% uncertainty in the calculated velocities. Similarily, the sta- 
tic pressure readings were examined and found to have an uncertainty of 
approximately 2% near the aft rake. This uncertainty in static pressure 
resulted in approximately a 5% uncertainty in the velocities at the aft 
rake. 
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The possibility that the total temperature profile at the aft rake 
was not that estimated in Figure 13 was also examined. A somewhat dif- 
ferent, but plausible, total temperature profile was plotted and is shown 
in Fig. 26 for 90% power. This change in the Ty profile resulted in a 
4.7% increase in the bulk static temperature and a 1.3% increase in the 
bulk static temperature and a 1.3% increase in the bulk velocity. These 
changes would result in a 4.8% decrease in the test cell mass flow rate 
(assuming engine flow rates were known). 

Finally, the assumption that the specific heat ratio was not 1.4 at 
the aft rake was examined. Changing y to 1.35 resulted in insignificant 
changes in the cell flow rate. 

The experimental data for the forward rake were also examined to dis- 
cover how uncertainties in static pressure and stagnation temperature 
profiles would affect the predicted velocities. Possible errors in the 
extrapolated pressure at the forward end of the pitot tubes resulted in 
an uncertainty of only 4.5% (Fig. 27) in the velocities. The velocity 
profile is sensitive to the assumed stagnation profile in much the same 
manner as discussed above for the aft rake. 

C. COMPUTER MODEL PREDICTIONS 

The pressure and velocity profiles that resulted from computer 
modelling are presented with the actual experimental data for each 
power setting in Figs. 10 through 12, 15 through 17, and 19 through 21. 

Sensitivity of the computer model predictions to specified test 
cell mass flow rate and engine mass flow rate were examined. A T0% de- 
crease in engine mass flow rate was found to alter the computer model 
results an insignificant amount. This is attributable to the lesser 
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contribution of the engine to the total flow rate through the augmentor. 

A 10% decrease in the specified test cell mass flow rate decreased the 
predicted velocity at the front rake a maximum of 10%. The maximum de- 
crease occurred near the augmentor wall, where the cell flow rate had the 
greatest influence. At the aft rake the 10% decrease in specified cell 
flow rate caused a maximum decrease in velocity of 7%, again near the 
augmentor wall. This decrease in specified flow rate caused the static 
pressure drop to decrease approximately 20%. These results are shown for 
the 90% power case in Figs. 28 and 29. 

The theoretical velocity profiles at the forward rake appear to be in 
good agreement with the experimental results. However, the experimental 
results showed less variations near the centerline than predicted. The 
predictions underestimated the centerline velocities at idle and 90% 
power and overestimated it at 100% power. 

The experimentally obtained aft rake velocities at idle power had 
considerable circumferential variations. The theoretical predictions, 
which are based on axisymmetric flow, appeared to be a reasonable esti- 
mation of the average behavior. In both the 90% and maximum power 
cases the computer model predicted a good velocity profile although the 
velocities were considerably greater (up to 63%) than the experimental 
values. This occurred to a greater extent at maximum power, indicating 
a possible decrease in predictive accuracy as mass flow rate was in- 
creased. Inaccurate specification of engine flow rates can significantly 
alter the cell flow rate input to the computer model. This in turn can 
significantly affect the predicted velocity profiles as discussed above. 
Inaccurate stagnation profiles can also contribute to errors in model 
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input data as discussed above. However, these possible inaccuracies 
cannot account for the entire difference in results obtained at full 
power. 

The theoretical pressure profiles (Figs. 10, 11, and 12) in all 
three cases indicate a large pressure drop at the augmentor tube inlet 
plane. The pressure decrease was greater at the augmentor wall than 
along the engine centerline . This theoretical drop in pressure was 
greater than the experimental pressure drop, the difference increasing 
with an increase in mass flow rate. The pressure rise from the minimum 
was in better agreement with the data but also was less accurate at 
higher power settings. The predicted behavior being much greater than 
experimental results is similar to the results of Walters and Netzer 
[Ref. 6*j using the computer code and to the results of Stevenson and 
Netzer [Ref. 1~] using the primitive variable code. In these earlier 
studies which compared predictions to experimental data in a one eighth 
scale test cell, good agreement in pressure profiles was obtained for 
zero engi ne-augmentor spacing. The discrepancy increased with increasing 
engine-augmentor spacing. One possible cause of this disagreement is the 
K-e turbulence model. It is being used in this application for the shear 
layer mixing of a practically unconfined jet surrounded by a very low 
velocity flow, both of which must then enter a confined diameter. 

Sensitivity of the computer model predictions to the magnitude of the 
effective viscosity was examined in a cursory manner by both increasing 
and decreasing the calculated values by a factor of 1.5. The computer 
model predictions were found to be fairly insensitive to these changes in 
calculated effective viscosity. For the case of a 50% increase in viscosity 
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Figs. 30 and 31 show less than 5% change in velocity at the front rake, 
15% at the aft rake, and no greater than a 6% decrease in static pressure 
drop along the augmentor wall. A 33% decrease in viscosity changed velo- 
city at the front rake less than 5%, 10% at the aft rake and decreased 
the static pressure drop along the augmentor wall no more than 9% (Figs. 
32 and 33). It appears that merely changing the magnitude of the effec- 
tive viscosity will not account for the differences in the predicted and 
measured pressure profiles. Further investigations of the two-parameter 
effective viscosity model are required. 

The sensitivity of the static pressure profiles to possible inaccura- 
cies in model inputs can account for much of the discrepancies at idle 
power but not for 90% and full power. 

Thus, it appears from both the velocity and pressure profile data 
that the predictive accuracy of the model decreases with increasing 
engine exit Mach number. 

Examination of the data from the velocimeters positioned above the 
TF-34 showed a test cell air velocity profile which increased from 
ceiling to engine. However, it has been shown by Speakman, Hayes, and 
Netzer jRef. 5"] that this condition has an insignificant effect on the 
predicted velocity and pressure profiles in the augmentor. 

The computer model appears to have the same capabilities and weak- 
nesses in predicting the flow in full scale test cells as it did for 
the subscale test cell. The subscale test cell, by design, provided 
nearly axi symmetric flow in the augmentor at all power settings and 
therefore somewhat better agreement was obtained for predicted velocity 
profiles far down the augmentor tube. Engine flow rates were also more 
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accurately known in the subscale tests. The most severe weakness of the 
model appears to be in the inability to correctly predict accurate sta- 
tic pressure profiles in regions where an expanding high velocity jet is 
captured by a fixed diameter tube. Velocity profiles and mixing rates 
appear to be reasonably predicted. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The two-dimensional, primitive-variable computer model proved to be 
quite good in predicting velocity profiles and mixing rates in the aug- 
mentor tube of a full scale jet engine test cell. The major weakness in 
this area was that the magnitude of the velocities far downstream were 
somewhat higher than actual with the discrepancy increasing with in- 
creasing engine exhaust Mach number (i.e. power setting). 

The greatest drawback to the computer model was its inability to 
accurately predict static pressure profiles in the augmentor entrance 
region where an expanding high velocity jet is captured by a fixed dia- 
meter tube. This inaccuracy also was found to increase with engine 
power level . 

If time and cost pose less restraints to future investigators it is 
recommended that the data gathering apparatus be expanded to include 
pitot-static and temperature measuring equipment at all locations where 
velocity data are desired. This would preclude the necessity for making 
the assumptions required in this investigation for stagnation tempera- 
ture profiles. 
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TF-34 Test Configurations 




Fig. 2 Pitot Rake Front View 
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Fig. 4 Pitot Rake Contruction 
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AUGMENTER TUBE 




Fig. 5 Rake Attachment 
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Fig. 6 Instrumentation 
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Fig. 7 Velocimeter Calibration Curves 
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Fig- 8 Pitot Tube/Thermocouple Labels 




Fig- 9 Static Port Labels 
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AUGMENTOR POSITION 



E 




I 



I 



I 



1 



l 



F 



/ 



/ 



o 



± 

o 



CO OJ CO 

0 — — 

1 « » 

£ 0I x^/N *dv 



J 

o 

CVJ 



50 



Fig. 10 Pressure Profile, Idle Power 



AUG mentor 




AUGMENTOR POSITION 



oo 



E 




in 



ro 



c\j — 




i 




i 



i 



i 



i 



i 



i 




o to o o m o ^ 

i 7" i cvj <nj ro 

• i i i 

£.°l X 2 UJ / N iiv 



Fig. 12 Pressure Profile, 100% Power 
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Fig. 13 Stagnation Temperature Profiles, Aft Rake 
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AXIAL VELOCITY , m/sec 
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AXIAL VELOCITY , m/sec 
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Fig. 17 Velocity Profile, Aft Rake, 100% Power 
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Fig. 18 Stagnation Temperature Profiles, Forward Rake 
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Fig. 21 Velocity Profile, Forward Rake, 100% Power 
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Fig. 22 Turbojet Test Cell Geometry for Computer Model 




Fig. 23 TF-34 Exhaust 
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Fig. 25 Average Velocity Profiles, Aft Rake 
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Fig. 26 Modified Stagnation Temperature Profile, Aft Rake, 90% Power 
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AU6MENT0R POSITION 
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Fig. 30 Sensitivity of Predicted Velocity Profile to 
Increase in Effective Viscosity, 90% Power 
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Fig. 32 Sensitivity of Predicted Velocity Profile to 
Decrease in Effective Viscosity, 90% Power 
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